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Similar to static systems, periodically driven systems can host a variety of topologically non¬ 
trivial phases. Unlike the case of static Hamiltonians, the topological indices of bulk Floquet bands 
may fail to describe the presence and robustness of edge states, prompting the search for new 
invariants. We develop a unified description of topological phases and their invariants in driven 
systems, by using scattering theory. We show that scattering matrix invariants correctly describe 
the topological phase, even when all bulk Floquet bands are trivial. Additionally, we use scattering 
theory to introduce and analyze new periodically driven phases, such as weak topologieal Floquet 
insulators, for which invariants were previously unknown. We highlight some of their similarities 
with static systems, including robustness to disorder, as well as some of the features unique to driven 
systems, showing that the weak phase may be destroyed by breaking translational symmetry not in 
space, but in time. 


I. INTRODUCTION 

Topologically non-trivial phases are characterized by 
the simultaneous presence of an insulating bulk and of ro¬ 
bust, conducting edge states.^Il^ The boundaries are pro¬ 
tected against localization by the system’s bulk proper¬ 
ties: its symmetries and the presence of a mobility gap. 
This link is termed bulk-boundary correspondence and 
it enables defining topological invariants, Z or Z2 val¬ 
ued quantities computed from the bulk system, which 
determine the presence and number of protected gapless 
modes at the boundary. 

The number and nature of topological invariants 
has been extensively studied for static Hamiltoni¬ 
ans. For non-interacting systems, there exist classifica¬ 
tions of symmetry protected topological insulators and 
super conduct or Their boundary states can be pro¬ 

tected not only by the three fundamental symmetries of 
the system: time-reversal, particle-hole, and chiral sym¬ 
metry, but also by symmetries of the underlying lat¬ 
tice. The latter are known as weak topological insula¬ 
tors in the case of translation symmetry or topological 
crystalline insulators for point group symmetries (reflec¬ 
tion, rotation, etc.), and can exist both in two and three 
space dimensions.Recent experiments have shown 
evidence for the existence and robustness of both.C^H^ 

Beyond static Hamiltonians, the topological phases of 
periodically driven systems have recently become the fo¬ 
cus of a wide range of researchPart of the appeal 
of so-called Floquet topological insulators is the possi¬ 
bility to ch ange a system’s properties by altering the 
driving field.While the physical properties of time- 
independent systems are mostly fixed during the fabri¬ 
cation process, a Floquet system can in principle host 
a variety of topological phases as a function of an ex¬ 
ternal, time-periodic perturbation. Examples include 
phases hosting Floquet Maj orana fermionsas well 
as robust chiral^^M^ helicalor counter-propagating 
edge modesPossible extensions of thes e idea s to in¬ 
teracting systems are also actively discussed. 

In spite of this intense activity, the topological na¬ 


ture of periodically driven systems is much less under¬ 
stood as compared to their time-independent counter¬ 
parts. Unlike wavefunctions of static systems, Floquet 
states are obtained from the unitary time-evolution op¬ 
erator over one driving period, the so-called Floquet op¬ 
erator, T. Each of its eigenstates accumulates a phase 
factor during one period T of the time evolution, as 
= exp(—i5T/fi)|'0), with 5 referred to as quasi¬ 
energy. While the Eloquet states and their associated 
quasi-energies are in many ways analogous to the eigen¬ 
states and energies of static systems, the periodic na¬ 
ture of the driving means that 5 is only defined modulo 
27rfi/T. The periodicity of the resulting Brillouin zone 
(BZ) in quasi-energy can lead to situations in which the 
original topological invariants fail to correctly capture 
the system’s behavior at an edge. Indeed, in Ref. |33| it 
was shown that robust chiral edge states can form in a 
system where all bulk bands are trivial, since the Eloquet 
operator in the bulk equals the unit matrix, IF = 1. 

Several works have made progress towar ds a full clas¬ 
sification of Eloquet topological insulatorsIn some 
cases, novel invariants have been formulated, which take 
into account the system’s driven nature. Eor strong two- 
dimensional (2d) Eloquet topological insulators, invari¬ 
ants have been found, bot h in t he presence and absence 
of time-reversal symmetry.l^^^^ Weak topological effects, 
however, have remained largely unexplored. Among 
the few examples to date, 2d chiral-symmetric systems 
hosting anomalous counter-propagating edge modes have 
been reported in Refs. l38ll^ and|42j In these works the 
existence and robustness to disorder of edge modes was 
identified, but the topological index responsible for their 
presence has remained up to now unknown. 

To accurately predict the possible non-trivial phases, 
the properties of Eloquet topological insulators - both 
weak and strong - need to be treated within a unified 
framework. To this end, rather than focusing on their 
difference with respect to static systems, we take ad¬ 
vantage of their similarities instead. In both types of 
phases, topological properties manifest themselves at the 
(quasi-)Eermi level, in the form of robust conducting edge 
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modes, which exhibit similar forms of spectral flow. As 
such, their non-trivial phases can be characterized in a 
unified manner, by using scattering theory. 

Previous works have considered the quantum transport 
properties of Floquet systemsxhe scattering prob¬ 
lem is typically defined in terms of a driven sample con¬ 
nected to multiple metallic electrodes and the scattering 
matrix relates the amplitudes of asymptotic scattering 
states in the metallic leads. In contrast, our main focus 
is not studying transport properties, but establishing a 
topological classification of Floquet systems in terms of 
the scattering matrix. To this end, rather than attach¬ 
ing metallic electrodes, we impose absorbing boundary 
conditions on a given sample to define a simplified, yet 
fictitious scattering problem.l^ As we will show, the re¬ 
sulting scattering matrix enables to fully determine the 
topological properties of periodically driven systems. 

The Fermi level scattering matrix has been used to for¬ 
mulate both strong and weak topol ogical invariants for 
time-independent Hamiltonians^^*^ as well as Id quan¬ 
tum walks,!^ having the advantage of being naturally 
tailored to the study of disordered systems. In this work, 
we extend this approach to periodically driven systems, 
showing different examples of its application. Specifi¬ 
cally, we formulate novel topological invariants for the 
models of Refs. |38] and [42l based on the constraints im¬ 
posed on the scattering matrix by the symmetries of the 
system. Furthermore, we show that the original scat¬ 
tering matrix invariants developed to characterize strong 
and weak static topological insulators can be readily ap¬ 
plied to driven systems, even in cases in which all bulk 
Floquet bands are trivial. This unified approach miti¬ 
gates the need for new topological indices, enabling the 
study of a wide class of both static and driven phases 
using the same invariant expressions. As an example, 
we turn to the model studied in Ref. [45] and reveal a 
richer topological structure consisting of both strong and 
weak phases. We showcase some features of weak phases 
which are unique to periodically driven systems, such as 
the possibility of gapping out the edge modes by breaking 
translation symmetry in time. 

The rest of our work is organized as follows. In Sec¬ 
tion [n] we study the chiral-symmetric Floquet system of 
Refs. |38| and |42l We determine the scattering matrix 
invariant and, based on its structure, predict the local¬ 
ization behavior of edge states in the presence of disorder. 
We find non-trivial phases of a Z 2 (as opposed to Z) na¬ 
ture, and conclude that the edge states are not robust to 
disorder which does not preserve chiral symmetry, con¬ 
trary to previous reports. We confirm this expectation 
by performing numerical simulations. In Section |m] we 
turn to the particle-hole symmetric Floquet topological 
insulator of Ref. EH which hosts both strong and weak 
topological phases. We find that the original, strong and 
weak scattering matrix invariants accurately describe the 
edge state robustness, even though all bulk bands are 
trivial. We conclude in Section lEl 


II. CHIRAL SYMMETRIC DRIVEN SYSTEM 

To date, there are few studies of weak topological ef¬ 
fects in periodically driven systems. Refs. [38] and |42| have 
observed this behavior in 2d systems possessing chiral 
symmetry. In this Section, we turn to one such model, 
the kicked quantum Hall system, and use scattering the¬ 
ory to determine its topological invariant. 

The model describes spinless fermions on a square lat¬ 
tice (lattice constant a = 1), where the position of each 
site is given by the vector r = UxX -h with Ux^y 
integers and x, y unit vectors pointing in the x- and y- 
directions, respectively. 

The tight binding Hamiltonian 

-f^KQH ^ ^ \x^x T {nx-)ny\-\- 

(l) 

—^C ^ \xix')'^y T T h.C., 

is expressed in terms of states \nx^ny) on lattice sites 
indexed by {nx^Uy). Here, Jx^y are the nearest neigh¬ 
bor hopping amplitudes in the x- and ^-directions. The 
system is placed in a uniform magnetic field, with a mod¬ 
eling the flux threaded in each plaquette. 

In the static case, Eq. § is the well-studied Hofstadter 
model,!^ showing a fractal pattern of gapped, Chern- 
insulating phases. As in Refs. [38] and HU we take Jx^y 
to be periodic functions of time, with Jx{t) = Jx and 
Jy{t) = Jy"^^S{t — m), m G Z. In other words, the 
hopping in the x-direction is kept constant, while the 
coupling in the y direction is turned on periodically with 
a period T = 1, but only at discrete times t = m. 

When the system is infinite in one or more direc¬ 
tions, or alternatively in the presence of periodic bound¬ 
ary conditions (PBC), momentum becomes a good quan¬ 
tum number, enabling us to express Eq. § in reciprocal 
space. The time-dependent momentum-space Hamilto¬ 
nian becomes 

"^KQH — Jx cos(A^3^) T Jy cos(^ky TixOL^ ^ ^ d{t m), ( 2 ) 

m 

where Ux labels sites in the magnetic unit cell. 

As long as there are no on-site potentials, the model 
shows chiral symmetry, expressed as a staggered gauge 
transformation. The chiral symmetry operator changes 
the sign of the wavefunction on one of the two sublattices 

r : \n^,ny) (-!)”-+”»|n^, Uy). (3) 

It is a unitary operator, = F^ = 1. In real space 
it anti-commutes with the Hamiltonian, FiZKQHF = 
—^KQH, while in reciprocal space it amounts to a tt shift 
in momentum, 

rHKQH(^tc, ky)^ = HKQH(^tc + TT, + Tt) 

= -^KQh(^cc, ^y)‘ 
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We set h = 1 and define the Floquet operator as prop¬ 
agating states from t = 0.5tot = 1.5: 


‘Fkqh = exp 


^1.5 

* / ' 
^0.5 


'Hkqh(^) dt 


(5) 


In this so-called symmetric time framep^^^ Eq. (§ can 
be written as 


‘Fkqh = (6) 


Owing to the chiral symmetry of the periodically mod¬ 
ulated Hamiltonian (§, the Floquet operator (§ also 
shows chiral symmetry, expressed as 


Tjr = exp(ma,7r) exp(in^7r). (7) 

While on the Hamiltonian level chiral symmetry amounts 
to changing the sign of the eigenvalues, in Floquet lan¬ 
guage it reverses the direction of the time-evolution 

rjT'^KQHFjT = ‘F;kQH' (^) 

We set a = 27r/3, such that the system has three bulk 
Floquet bands. In the following we consider Eq. ® ei¬ 
ther in a strip geometry (infinite in finite in x), or 
discretized on a square lattice of L x IF sites. All numer¬ 
ical results are obtain ed for tight-binding models defined 
using the Kwant code.^^^^ 

Typical banstructures obtained by diagonalizing the 
Floquet operator JicqhI'^) = exp{—ie)\^) in an infinite 
strip geometry are shown in Fig. As a function of the 
parameters Jx^y, the model exhibits a variety of topolog¬ 
ical phases. The top and bottom quasi-energy gaps of 
Fig. 0 {s = 7t/2 and 5 = 37r/2) show one or more edge 
modes with a non-zero net chirality, a feature reminiscent 
of the quantum Hall effect. The boundary modes occur 
on all edges, irrespective of orientation, and have been 
understood^ in terms of the Floquet winding number of 
Ref. [331 

We focus on edge modes in the gap at 5 = tt which 
have no net chirality, but come in counter-propagating 
pairs. Each pair is separated by a momentum difference 
of TT. Furthermore, unlike boundary states in the top and 
bottom gaps, the counter-propagating modes depend on 
the orientation of the edge, similar to weak topological 
insulators. For the values of Jx^y used in the left panel of 
Fig. there are no gapless edge states on boundaries par¬ 
allel to either the x- or ^-direction. In the middle panel 
both types of boundary show one pair of edge modes, 
while in the right panel there are two pairs on the edge 
parallel to and none along x. There is to date no 
known topological invariant responsible for the presence 
and robustness of these edge states. 

The presence of counter-propagating modes has been 
previously identified as being a consequence of chiral 
symmetry. Indeed, in a strip geometry the relation © 
becomes 

^jrJ~¥.CiR{ky)Tj: = ~ (^) 


such that for each state at energy 5 and momentum ky 
there must be a state at —5 and ky — it. Therefore, edge 
modes at e = Ttt necessarily come in pairs. In the follow¬ 
ing we will determine the nature of this topological phase 
as well as the explicit form of its topological invariant. 

In the absence of disorder momentum is a good quan¬ 
tum number and a single pair of edge states is protected, 
since the states must be separated by a momentum dif¬ 
ference Aky = TT. Figure also shows a situation in 
which there are two pairs oicounter-propagating modes 
at £ = TT, raising the question of whether this is a dis¬ 
tinct phase or whether different pairs of edge modes can 
annihilate and gap out if they overlap. The answer will 
determine whether the phase is of a Z type, with any 
number of protected edge modes, or of a Z 2 type, where 
only an odd number of pairs is protected. We test this 
hypothesis in a controllable fashion, by considering two 
copies of Eq. 0, coupled in such a way as to preserve the 
chiral symmetry 0 of the combined system (see Fig.[^. 
Each of the two subsystems has one pair of edge modes, 
but one of them is shifted in momentum relative to the 
other, ky ^ ky 6. This enables us to effectively slide 
one pair of edges in the quasi-BZ by changing the value of 
(5. In a strip geometry, the Hamiltonian of the combined 
system has the form 


-f^2xKQH = F ( V F 1”’®’-^' + 


i=l,2 


J, 


^ inxO!.-\-ik. 




^ +h.c.y 

i/i' / 


( 10 ) 


where the index = 1,2 labels the two copies, and Jc is 
the time-independent coupling between them. In Fig. 
we show that as the two pairs of edge modes are made to 
overlap, gaps open in the spectrum. Therefore, there are 
only two topologically distinct edge mode configurations, 
one pair of counter-propagating modes versus no edge 
modes, indicating that the non-trivial phase is protected 
by a Z 2 , as opposed to a Z invariant. 

To construct this invariant as well as test the conduct¬ 
ing properties of the edge, we turn to the scattering ma¬ 
trix formalism. For a static Hamiltonian discretized on 
a finite lattice, transport properties are usually obtained 
by connecting one or more infinite, translationally in¬ 
variant leads to the system. This enables to compute the 
scattering matrix, which in the two-lead case reads 

«=(?(,). (11) 


in terms of the reflection and transmission amplitudes of 
lead modes, and 

For periodically driven systems, the scattering matrix 
associated to a Floquet operator can be computed in 
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FIG. 1. Bandstmcture of the Floquet operator ^ in an infinite strip geometry (infinite in y, L = 21 sites in the x-direction). 
Parameters are a = 27tI3, Jx = 27r/3, and Jy = tt, 27r,37r, in the left, middle, and right panels, respectively. Red (dark gray) 
indicates states localized on the left boundary and blue (light gray) states on the right. From left to right, the bulk gap at e = tt 
shows zero, one, or two pairs of counter propagating edge modes, separated by a momentum difference Aky = tt (horizontal 
arrows). 



FIG. 2. Side view of the Hamiltonian (10), consisting of two 
copies of 0 labeled by j = 1, 2 (top and bottom layers). One 
of the three unit cells shown is marked by a box. The coupling 
within the top and bottom layers are marked with orange 
and red lines, respectively, while the one between the layers 
is shown in blue. Notice that chiral symmetry is preserved, 
as there is no coupling between layers in the same unit cell. 


a slightly different fashion.On a square lattice of 
L X W sites, the Floquet operator (|^ is a LW x LW 
matrix. Rather than attaching infinite leads, we intro¬ 
duce absorbing terminals at the two sides of the system, 
= 1 and rix = L. Defining a 2W x LW projector onto 
these terminals. 


P = 


1 if G {1,1/} 
0 otherwise 


( 12 ) 


we can express the quasi-energy dependent Floquet scat¬ 
tering matrix through the formula 


S{e) = -P [1 - ^(l - P^P)] ^ e'^PP"^, 


(13) 


where the superscript T stands for transposition. The ex¬ 
pression (13) can be understood by expanding the inverse 
matrix in a geometric series, in which each subsequent 
term describes time-evolution over an additional period. 
Time evolution stops after the state reaches the absorb¬ 
ing terminals, P, and continues otherwise, 1 — p^ppH 
The absorbing terminals only act stroboscopically, at the 
beginning and end of each time period, such that Eq. (13) 


describes a simplified, fictitious scattering problem. Nev¬ 
ertheless, the unitarity of the Floquet operator P implies 
that the scattering matrix is also unitary, S{s)S\e) = 1, 
which can be verified though a direct calculation. 

In this setup, the Floquet scattering matrix takes the 
form of Eq. 0), enabling to compute the transmission 
from one side of the system to the other, G = LitH. At 
5 = TT, the presence of one pair of counter-propagating 
edge states leads to a quantized transmission G = 2, since 
there is one right- and one left-mover on each edge. 

Eirst we identify the constraints imposed by the chiral 
symmetry on S. To this end, we consider a system with 
an odd number of sites in the y direction {W odd) and ap¬ 
ply twisted boundary conditions, as W) = 0). 

The Eloquet operator, and therefore also the scattering 
matrix, now become functions of the twist angle 0, which 
plays the same role as momentum in the chiral symme¬ 
try constraint of Eq. Plugging this constraint into 
the Eloquet scattering matrix definition (13) leads to a 
relation 


rs5(£,<^)r5 = 5t(-e,<^-7r), (14) 

where F 5 'isa 21 Ex 21 E matrix defined by F^ = PTjrP^. 
In other words, F^ changes the sign of one sublattice in 
the absorbing terminals at = 1,P. 

The relation 0 is reminiscent of that found for chiral 
symmetric static systems ,1^ where at zero energy there 
exists a basis in which the scattering matrix is hermi- 
tian. For the periodically driven system studied here, 
the TT momentum shift induced on the Floquet level car¬ 
ries over in the scattering matrix description. Therefore, 
we introduce the chiral basis S = E^P, in which the re¬ 
flection sub-block obeys r((/)) = (0 — tt) at e = tt. This 

enables us to formulate a topological invariant, by noting 
that for this quasi-energy 

det r(0) = det* r((/) — tt). (15) 
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FIG. 3. Bandstmctures of the model (10) in the same geometry as Fig. We set a — Jx — 27r/3, Jy — 27r, and Jc = 0.6. 
In the left, middle, and right panels, 6 = O.Ttt, O.SStt, and O.OStt, respectively. Sliding one pair of edge modes on top of the 
other (horizontal arrows) causes them to annihilate and opens gaps in the spectrum. This indicates they are protected by a Z2 
invariant. 



FIG. 4. Two possible ways in which the phase of det r can 
evolve as the twist angle is advanced from 0 = 0 to 0 = 
TT. The solid and dashed lines indicate topologically different 
scenarios, since they cannot be deformed into each other as 
long as the bulk mobility gap remains open (det r 7^ 0). The 
presence of protected 7r-crossings signals a topologically non¬ 
trivial phase in our choice of gauge. 


By Eq. (15), the complex phase of the determinant of 
r must change sign as the twist angle 0 is advanced by tt. 
Since the complex phase must be a continuous function 
of 0, only two scenarios are possible. In the interval 
4> G [0,7r] (half of the effective BZ), either the phase of 
det r crosses 0 an odd number of times, or it crosses tt 
an odd number of times, as shown in Fig. 

The odd number of 0-crossings or 7r-crossings in half of 
the BZ is a Z2 topological invariant, and cannot change 
as long as there is a mobility gap in the bulk: det r 7^ 0. 
While the two scenarios are topologically different, which 
of them occurs is a matter of definition. In order to pre¬ 
serve the form of the chiral symmetry requirement (|^, 
we have used an odd number of sites in the ^-direction, 
making r a matrix with odd dimensions. Therefore, de¬ 
pending on which of the two sub-lattices changes sign un¬ 
der the application of F^, the determinant of r will also 
change sign, turning 0-crossings into 7r-crossings, and vice 
versa. For definiteness, we chose a gauge in which the chi¬ 
ral basis has det F^ = — 1 and identify 7r-crossings with a 


topologically non-trivial phase. Choosing the other sign 
would imply interchanging the labels of Fig. 


By numerically computing the reflection matrix as a 
function of twist angle, we have checked that this Z2 in¬ 
variant correctly describes the topological phases of the 
model. Furthermore, the relation © still holds when 
disorder is added to the system, provided the latter does 
not break chiral symmetry. We therefore expect the edge 
state transmission to be robust upon the inclusion of ran¬ 
dom hopping strengths Jx,y On the other hand, the 
system should fully localize if we add random on-site po¬ 
tentials, which violate the chiral symmetry requirement 
of Eq. and ( p!^ . 


We confirm this numerically by including disorder in 
one of two different ways. For random bond disorder, 
we substitute Jx,y ^ with 6j drawn in¬ 

dependently for each bond from the uniform distribu¬ 
tion [—t/, t/], with U the strength of disorder. In the 
case of on-site disorder, we introduce a random, time- 
independent chemical potential term to the Hamiltonian 
0, where is drawn independently 

for each lattice site from the same, uniform distribution. 
The transmission scaling results of Fig. confirm our 
expectations. Bond disorder leads to an algebraic de¬ 
cay of the edge transmission, G ^ characteris¬ 

tic of two-dimensional statistical topological insulators .1^ 
The edge remains deloc alized, being pinned to a one¬ 
dimensional critical point and the phase of det r 

still shows protected 7r-crossings. In contrast, we find 
that on-site disorder leads to an exponential suppression 
of edge transmission, G ^ exp(—cL), c = const., signal¬ 
ing localization. 


Before concluding this Section, we analyze the top and 
bottom energy gaps appearing in Fig. Away from 
5 = 0, TT, the system does not obey chiral symmetry, and 
there is an imbalance in the number of left and right 
movers on an edge. Since these edge states manifest in 
the absence of any symmetries, the scattering matrix in¬ 
variant originally developed to describe the Chern num- 
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FIG. 5. Log-log plot of the edge transmission of model § at 
£ = TT as a function of system length L for bond (red) and 
on-site (blue) disorder, with W = 40, a = Jx = 27r/3, and 
Jy — 27 t. For on-site disorder we have used a disorder strength 
U =1.8, while for bond disorder it was set to = 0.45. In 
each case 600 — 800 independent disorder realizations were 
averaged over. The solid and dashed lines show the expected 
algebraic and exponential decay of transmission, respectively. 



L 


FIG. 6. Setup and driving protocol of Hamiltonian |T^ . We 
consider an L x IF hexagonal lattice (sublattices A and H), 
with Bravais vectors ax and ay . When Jn = 0 only one of the 
hopping amplitudes Jx,y,z (blue, green, and red, respectively) 
is turned on at any given time. For Js — 2 particles at any 

site in the bulk return to their original position after two driv¬ 
ing periods, while those on the edges form chiral propagating 
modes, as shown by the arrows. 


ber of static system^^^^ correctly captures the net chi¬ 
rality of the edge states. By applying twisted boundary 
conditions as before, we write the topological invariant 
as the winding number of det r 


CSM = ^. dct> — ] 
Zttz Jq dcp 


For the parameters in the left and right panels of Fig. 
we find that the top gap has a scattering matrix Chern 
number Csm = +I 7 while the bottom one has Csm = ~ 1 - 
By repeating the calculation for the parameters of the 
middle panel, we find Csm = — 2 and Csm = +2 for the 
top and bottom gaps, respectively. 

As we have shown, scattering theory enables to fully 
characterize the topological phases of the model The 
counter-propagating modes appearing in the gap at e = tt 
can be understood in terms of an invariant readily ob¬ 
tained by means of a symmetry analysis. Chiral modes 
present in other gaps are accounted for by the invariant 
(16), originally developed in the context of static Chern 


insulators. In the following Section we will further ex¬ 
plore the connection between the scattering matrix in¬ 
variants of static systems and the non-trivial phases of 
Floquet topological insulators. 


III. PARTICLE-HOLE SYMMETRIC DRIVEN 
SYSTEM 

One of the advantages of scattering theory in the study 
of periodically driven topological phases is that it allows 
for the study of both time-independent and Floquet topo¬ 
logical insulators within the same, unified framework. In 
this Section we use the scattering matrix formalism to 


reveal and characterize the rich structure of a model con¬ 
sisting of both strong and weak topological phases. As 
we will show, scattering matrix invariants of Floquet sys¬ 
tems can take the same form as in their static counter¬ 
parts, even when bulk Floquet bands fail to capture the 
non-trivial character of a phase. 

We consider an example of particle-hole symmetric 
driven system, introduced in Ref. [45l It describes spin¬ 
less fermions on a hexagonal lattice in the presence of 
nearest neighbor hopping. The unit cell consists of two 
sites belonging to the two sub-lattices, A and B, and the 
Bravais vectors of the lattice ax and ay are defined as in 
Fig. in 

Setting the lattice constant to 1 , the momentum space 
tight-binding Hamiltonian has the form 

^ (‘^A.kT.k) (cB,k) ’ 

k \ ’ / 

with 

n ={Jx cos{ky - kx) + Jy cos{ky) + Jz)<Jx- 

{Jx sm{ky - kx) + Jy sin{ky))ay. 


Here, ^ are fermionic creation and anni¬ 

hilation operators on the A and B sub-lattices and Jx,y,z 
are independent hopping amplitudes describing three dif¬ 
ferent types of bonds, as shown in Fig. The Pauli 
matrices appearing in B parametrize the sub-lattice 
degree of freedom. 

In the absence of driving, the Hamiltonian (18) belongs 
to class BDI in the Alt land-Zirnbauer classification.l^ It 
is characterized by time-reversal, particle-hole, as well 
as chiral symmetries, all of which square to +1. The 
corresponding operators are T = 1C, V = (TzIC, and C = 
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az, with K, complex conjugation, such that 

n{k) =H*{-k), 
aznik)az = -n*{-k), 
az'H{k)(Jz = - 'H{k). 


(19) 


We choose a driving protocol which involves the cyclic 
modulations of the hopping amplitudes J^P^Each of the 
three hoppings contains a uniform term which is in¬ 
dependent of time r, as well as a term Jg which is peri¬ 
odically varied in a stroboscopic manner. The three-step 
driving protocol reads: 

1- Jx — Js Ju^ ~ Ju 

for nT < r < nT + T/3, 

2. Jy = t/g “h Ju^ Jz,x ~ Ju 

for nT + T/3 < T < nT + 2T/3, 

3- Jz ~ Js T Ju-) Jx,y ~ Ju 

for nT + 2T/3 < r < nT + T, 

with T the driving period. 

In the following we set T = 3 and h=l expressing the 
Floquet operator as the product 

T = exp(—exp(—i?^ 2 ) exp(—(20) 

where 1-ii are the Hamiltonians during the three steps 
of the driving protocol shown above. Even though at ev¬ 
ery instance of time the Hamiltonian (18) is time-reversal 
symmetric, the sequence in which the x, and z-type 
hoppings are modulated implies that the Floquet opera¬ 
tor (20) has broken time-reversal symmetry. 


In the simplest case we set = 0, such that during 
each of the three steps of the driving protocol only one 
of the Jx,y,z hoppings is non-zero. When additionally 
Js = 7r/2 a particle is transferred with unit probabil¬ 
ity between neighboring sites. Therefore, particles in the 
bulk fully encircle one hexagonal plaquette in two driving 
periods, leading to the formation of dispersionless (flat) 
bulk bands. In the presence of a boundary, the same 
driving protocol leads to the formation of a chiral propa¬ 
gating mode on the edge of the lattice, as shown in Fig.|^ 
We recover the dispersionless bulk bands as well as the 
chiral edge modes in the quasi-energy banstructure of the 
system (see Fig. [7[i). 

The emergence of the flat bulk bands attracts signif¬ 
icant attention on its own right, in the context of pos¬ 
sible realizations of exotic many-body phases in such 
systemsIn the case of a non-trivial topology of the 
flat band, reflected in non-zero band Chern number, chi¬ 
ral edge states can emerge at the boundary of a finite 
system as a result of bulk-boundary correspondence. 

Interestingly, here for = 0 and Js = 7r/2 the bulk 
Floquet operator, T = i{cos{kx)<Jx + sin(/ca,)cr^) is inde¬ 
pendent of ky, meaning that both of the flat bands have 



TT -TT A; 0 
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FIG. 7. Bandstructures of the model ( F8| at Js = 7r/2 in 
a strip geometry (infinite along ax^ W — 20). In panel (a), 
Ju = 0 and there are two flat bulk bands at quasi-energies 
£ = ±7r/2 as well as one chiral edge mode present at all quasi¬ 
energies. In panel (b) we set Ju = 0.25, such that there is one 
chiral mode at £ = tt and a pair of counter-propagating modes 
at £ = 0. The strong and weak invariants at £ = tt read 


(Gsm, I 


y, 7 r) = (—1,1,1) in both panels. At e: = 0 the 


system in panel (a) is in a (—1,—1,—1) phase, while panel 
(b) has (0,—1,—1). The color scale denotes the eigenstate 
intensity on the first and last 10% of lattice sites. 


zero Chern number and all bulk states are trivially local¬ 
ized. The topological protection of the edge states can 
nevertheless be deduced from the topological invariant 
formulated in Ref. |33l which takes into account the full 


time evolution of the Hamiltonian (18) throughout the 
driving cycle. 

In the following we show that scattering theory can 
not only reproduce this invariant, but reveals a much 
richer topological structure. To this end, we consider 
finite systems of size LxW (see Fig.[^, where L denotes 
the number of vertical bonds in the ax direction and W 
the number of zig-zag chains in the ay direction. The 
scattering matrix associated to the Floquet operator is 
obtained by introducing absorbing terminals on the first 
and last zig-zag chain. 

We repeat the analysis of Section [TTj finding that the 
presence of a chiral edge mode leads to a quantized trans¬ 
mission G = 1. As before, we apply twisted boundary 
conditions (twist angle 0) in the ax direction and express 
the strong topological invariant as the winding number 
of det r, Eq. Remarkably, the scattering matrix in¬ 

variant originally developed to compute the Chern num¬ 
ber of static systems correctly describes this phase, even 
though all bulk bands are trivial. We find that the scat¬ 
tering matrix Chern number equals Csm = —1 foi* all 
quasi-energies 5 in the bulk gaps of Fig. 

Beyond the topological classification in terms of the 
strong invariant, the Hamiltonian ( p!8| ) shows a richer 
structure due to the presence of particle-hole symmetry. 
The cyclic nature of the driving protocol leads to a Flo¬ 
quet operator (20) which breaks time-reversal symmetry. 


but shows particle-hole symmetry (PHS) of the form 

^(k)=a,^*(-k)a,. (21) 

As such, for every eigenstate at quasi-energy 5 and 















momentum k there must be a state at —5 aud — k. A 
chiral edge mode preseut at 5 = 0, ±7r must therefore 
exist at poiuts where k = — k. This is the case for the 
baudstructure showu iu Fig. where the edge state 
crosses = 0 for e = tt aud kx = tt for s = 0. The two 
edge mode coufiguratious at 5 = 0 aud e = tt cauuot be 
coutiuuously deformed iuto each other without breakiug 
particle-hole symmetry or closiug the bulk gap, sigualiug 
the two phases are topologically distiuct. 

Iu other words, the parity of edge modes at momeuta 
kx = 0, TT appeariug iu Fig. is topologically protected. 
Iu a static topological supercouductor this protectiou 
would be expressed iu terms of a weak iuvariaut) ^^ * ^^ * ^^ * ^ 
Motivated by this fact, we aualyze the gapped phases of 
the model ( [2Q| ) iu terms of the scatteriug matrix iuvari- 
auts developed to classify weak topological supercouduc- 
tors. 

To ideutify how the iuvariaut emerges we aualyze the 
symmetries coustraiuiug the scatteriug matrix. Pluggiug 
the particle-hole symmetry coustraiut (21) iuto the scat¬ 


teriug matrix expressiou (13) leads to a relatiou 


S{s) = azS''{-e)az 


( 22 ) 


where az acts ou the sub-lattice degree of freedom iu 
the absorbiug termiuals. As a cousequeuce of Eq. 
at the particle-hole symmetric quasi-euergies 5 = 0,7r 
there exists a basis iu which the scatteriug matrix, aud 
therefore also its reflectiou sub-block are real. We choose 
the basis f = UrU\ with U = diag(l, i, 1, i, ..., 1, i) such 
that f = f*, aud iutroduce Z 2 weak scatteriug matrix 
iuvariauts couutiug the parity of edge modes ai kx = 0 
aud TT, 


= sign detf(()) = 0 ), 
i^x, 7 r = sign detf((() = tt). 

Here, 0 = 0, tt correspouds to applyiug periodic or auti- 
periodic bouudary couditious iu the ax directiou, respec¬ 
tively. The weak iuvariauts iu the ay directiou, 
cau be defiued iu a similar mauuer, by iutroduciug ab¬ 
sorbiug termiuals at = 1,1/ aud usiug (auti-)periodic 
bouudary couditious aloug ay^ 

The topological uumbers of Eq. ( [T^ aud ( [^ have 
features similar to those defiued for static systems. The 
stroug Z iudex couuts the uet uumber of chiral modes, 
whereas the weak Z 2 iuvariauts determiue the parity of 
modes crossiug at momeutum 0 or tt. As such, the weak 
iuvariauts are coustraiued by the parity of the stroug iu¬ 
dex CsM as = (—1)^^^, with i = Charac- 

teriziug a topological phase therefore requires specifyiug 
C'sM, as well as i^y^n^ siuce their values cau chauge 
iudepeudeutly.l^ Iu accord with the termiuology used iu 
static weak topological supercouductors, we label a phase 
where au odd uum ber of edge modes cross momeutum tt 
as uou-trivialj ^^ l ^^ l 

Eor Ju = 0 aud Jg = 7r/2 we hud that iu both di- 
rectious i = x^y we have = “1 aud = 1 at 


5 = TT, while their sigus are reversed at 5 = 0. There¬ 
fore, the 5 = 0 phase of Eig. is uou-trivial iu both 
a stroug aud a weak seuse. Eveu though the stroug iu¬ 
variaut CsM = — 1 remaius uuchauged, weak iudices cau 
distiuguish betweeu the two phases, such that au iuter- 
face betweeu the phase at 5 = 0 aud 5 = tt would host 
a pair of couuter-propagatiug edge modesA purely 
weak phase is obtaiued iu Eig. [^, where turuiug ou a 
small uuiform hoppiug, = 0.25, chauges the stroug 
iuvariaut to Csm = 0, while keepiug ^ 7 ^ = —! at 5 = 0. 

As iu time-iudepeudeut topological super couductors, 
the weak iuvariauts iutroduced above are protected by 
the combiuatiou of PHS aud trauslatioual symmetry of 
the system. Breakiug the latter, for iustauce by iutroduc¬ 
iug a staggered modulatiou of hoppiug amplitudes, would 
fold the quasi-BZ of Eig. [^iu momeutum, couvertiug the 
uou-trivial weak iuvariauts at 5 = 0 iuto trivial oues. 
Nevertheless, we hud that weak iuvariauts are robust to 
disorder as loug as the latter preserves trauslatioual sym¬ 
metry ou average, aud the 5 = 0 phase iu Eig. cau be 
thought of as the periodically driveu aualogue of a statis¬ 
tical topological iusulator ouce disorder is iutroduced.l^ 
Iu this phase, we verify that addiug a raudom compo- 
ueut to Ju^ drawu iudepeudeutly for each boud from a 
uuiform distributiou, leads to au algebraically decayiug 
edge trausmissiou G ^ The couuter-propagatiug 

edge states avoid localizatiou due to a combiuatiou of 
particle-hole symmetry aud average trauslatioual sym¬ 
metry, effectively formiug a Eloquet Kitaev edge.^ 

There are however importaut differeuces betweeu the 
topological classificatiou of time-iudepeudeut supercou- 
ductors aud that of the model ( [^ . Eirst, iu periodically 
driveu systems there are two PHS quasi-euergies at which 
the weak iuvariauts cau be defiued, 5 = 0 aud tt, allow- 
iug for a richer topological structure. Eveu though the 
weak iuvariauts at 5 = 0, tt are protected by the same set 
of symmetries, PHS aud trauslatiou aloug either ax or 
a^, they cau vary iudepeudeutly of each other. This is 
to be co utrast ed to the topological classificatiou of static 
systems ,l^^^^iu which differeut iuvariauts always require 
differeut sets of protectiug symmetries. 

Auother feature uuique to weak topological Eloquet 
systems cau be traced back to the periodicity of quasi- 
euergy. Iu static systems a weak phase may be destroyed 
by a dimerizatiou iuduced breakiug of trauslatioual sym¬ 
metry, leadiug to a doubliug of the uuit cell aud a foldiug 
of the BZ iu momeutum. Iu periodically driveu systems 
there is au extra directiou iu which the BZ cau be folded, 
that of quasi-euergy. Such a breakiug of time transla¬ 
tion symmetry cau be achieved by iutroduciug a half- 
frequeucy compoueut iu the driviug protocol, effectively 
doubliug the driviug period. This leads to a uew phase 
at 5 = 0, with iuvariauts (Csm, ^cc, 7 r 5 obtaiued by 
the compositiou of iudices iu the origiual (pre-foldiug) 
5 = 0 aud 5 = TT phases. While iu our example the folded 
phase is obtaiued simply by the Z additiou of the stroug 
iudex aud the Z 2 additiou of the weak oues, more com¬ 
plex compositiou rules may apply iu a geueric settiug.^^ 
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FIG. 8. Bandstmcture of the model (18) in an infinite strip 
geometry with boundaries along ax (top, W = 30) and ay 
(bottom, L = 30) directions. We use Jg = 7r/2 and in¬ 
equivalent uniform components of the hopping amplitudes, 
Ju,x = —0.05, Ju,y — Ju,z = 0.5. In the right panels 
the BZ is folded in quasi-energy by changing the strobo¬ 
scopic part of Jx to Jg^x = 7r/2 + 0.6 on every second pe¬ 
riod. The invariants of the pre-folding phases (left panels), 
(Gsm, J^x, 7 T,J^y, 7 T) = (0, -1, -1) at e: = 0 and (0, -1,1) at e = tt 
combine to form a new gapped phase with (0,1, —1) at £ = 0. 
The color scale denotes the eigenstate intensity on the first 
and last 10% of lattice sites. 


We show an example of time folding in Fig. ob¬ 
tained using inequivalent hopping amplitudes Jx,y,z in 
Eq. (18). Setting Jg = 7r/2 and the uniform components 
to Ju,x = —0.05 and Ju,y = Ju,z = 0.5, we obtain phases 
with indices (0, —1, —1) at 5 = 0 and (0, —1,1) at 5 = tt. 
The folding is implemented by repeating the three-step 
driving protocol described above twice, and altering the 
stroboscopic component of Jx during every second pe¬ 
riod, as Js^x = Js with 5Js = 0.6. This leads to 

a new gapped phase at 5 = 0, with invariants (0,1, — 1), 
obtained by the composition of the original ones. 


IV. CONCLUSION 


us to reveal a richer topological structure in a previously 
studied model (Section [nl|. We wrote down expressions 
for both strong and weak topological invariants, showing 
they correctly describe the non-trivial phases even when 
all bulk bands are trivial. 

When driven systems show weak phases for which no 
time-independent counterpart is known, as in Section |TI| 
scattering theory can be readily used to formulate novel 
invariants, based on the constraints imposed on the scat¬ 
tering matrix by the symmetries of the system. Since 
scattering matrix invariants are naturally tailored to the 
study of disordered systems, we were able to use their ex¬ 
pressions to deduce the conditions under which the weak 
phase survives disorder. 

Finally, we have shown that breaking translational 
symmetry in time can be used to alter the topological 
invariants of a phase. This novel feature is unique to 
Floquet topological insulators and allows a greater level 
of control in the design and manipulation of topological 
phases. Using the model of Section [ITT] we have demon¬ 
strated this idea and altered the topological invariants of 
a weak phase by doubling the driving period. 

We emphasize that, even though we have focused on 
two specific models, our scattering matrix approach to 
Floquet systems is completely general. It can, for in¬ 
stance, also be applied to time reversal invariant Floquet 
topological insulators or higher dimensional systems, as 
was the case for the scattering matrix invariants of static 
systemsOur method is numerical in nature and re¬ 
quires systems large enough to avoid finite-size effects 
and show a well-defined bulk mobility gap. However, 
given a system size, Floquet scattering matrix invariants 
are more efficient to compute than ones based on wave- 
functions. No knowledge of Floquet eigenstates or time- 
integration is required, and the biggest computational 
cost comes from inverting the time-evolution operator in 
Eq. 

We hope our work will pave the way towards a full 
classification of topological phases in periodically-driven 
systems, which would parallel that of time-independent 
Hamiltonians. Additionally, our work will prove useful 
in determining the robustness to disorder of non-trivial 
phases, an essential ingredient for their experimental re¬ 
alization. 


In many static systems, topological properties are fixed 
during the fabrication process. Part of the active inter¬ 
est in Floquet systems comes from the possibility to tune 
these properties by altering the driving protocol. Nev¬ 
ertheless, the topological classification of periodically- 
driven systems has remained largely unexplored. 

We have shown that topological phases of periodically- 
driven systems can be analyzed in a unified framework, 
by using scattering theory. In some cases, both static 
and Floquet systems may be characterized by the same 
scattering matrix invariants, owing to the similar trans¬ 
port properties of their edge modes. This has enabled 
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